On the corona of magnetars 

Yury Lyubarsky and David Eichler 

Physics Department, Ben-Gurion University of the Negev, Beer-Sheva, Israel 



ABSTRACT 

Slow dissipation of non-potential magnetic fields in the magnet osphere of the 
magnetar is assumed to accelerate particles to hundreds MeV along the mag- 
netic field lines. We consider interaction of fast particles with the surface of the 
magnetar. We argue that the collisionless dissipation does not work in the at- 
mosphere of the neutron star because the two-stream instability is stabilized by 
the inhomogeneity of the atmosphere. Rather, the dominant dissipation mech- 
anism is collisional Landau level excitations followed by pair production via the 
deexcitation gamma-rays ultimately leading to electrons with the energy below 
the Landau energy. We show that, because of the effects of the superstrong mag- 
netic field, these electrons could emerge from the surface carrying most of the 
original energy so that a hot corona arises with the temperature of 1 4- 2 MeV. 
This extended corona is better suited than a thin atmosphere to convert most of 
the primary beam energy to non-thermal radiation and, as we show, most of the 
coronal energy release is radiated away in the hard X-ray and the soft gamma- 
ray bands by Comptonization and bremsstrahlung. The radiation spectrum is 
a power-law with the photon index 1 < a < 2. The model may account for 
the persistent hard X-ray emission discovered recently from the soft gamma-ray 
repeaters and anomalous X-ray pulsars and predicts that the radiation spectrum 
is extended into the MeV band. 

Subject headings: instabilities - plasmas - stars:magnetic field - stars:neutron 



1. Introduction 

It is now commonly accepted that Soft Gamma Repeaters (SGRs) and Anomalous 
X-ray Pulsars (AXPs) are magnetars, neutron stars with extremely high magnetic fields, 
B ~ 10 14 -4 10 15 G (Duncan & Thompson 1992; Paczynski 1992; Thompson & Duncan 
1995). Activity of these sources is fed by the energy of this magnetic field (see Woods & 
Thompson 2004 for a recent review). It was found recently (Kuiper et al. 2004, 2006; 
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Mereghetti et al. 2005, Molkov et al. 2005) that the persistent pulsed X-ray emission from 
these sources has a nonthermal spectrum extending up to > 100 keV. The luminosity of this 
tail, L ~ 10 36 erg/s, exceeds the thermal luminosity from the star's surface. The spectra are 
exceptionally hard with photon indices typically in the range 1 -f- 2 so that the luminosity 
peaks above 100 keV. In this paper, we discuss origin of this emission. 

According to Thompson, Lyutikov & Kulkarni (2002), magnetic field in magnetar's mag- 
netosphere is generally non-potential. Slow untwisting of the magnetic field lines generates 
the electric field; particles are accelerated in this field until they produce electron-positron 
pairs filling the magnetosphere and providing charge carriers for the electric current (Be- 
loborodov & Thompson 2007). Pairs are formed from gamma-photons produced by cyclotron 
scattering of the surface radiation of the neutron star on primary particles. In magnetar's 
magnetic field, the energy of the first Landau level is high, 

e B = m e c 2 ^l + 2B/B g « 3.5 v^MeV; (1) 

so that the scattered photon is immediately converted into an electron-positron pair. Here 
B = 10 15 i?i5 G is the magnetic field of the star, B q = m 2 e c^ /he = 4.4 x 10 13 G. If the surface 
emits radiation with the characteristic photon energy e = lOeio keV, the cyclotron scattering 
occurs when the Lorentz factor of the electron reaches the value 

i= B_n h l = i ^ (2) 
B q e £io 

Note that the recoil effect should be taken into account in this estimate (Thompson & 
Beloborodov 2005). Eventually the energy is deposited in pairs with the Lorentz factor 
something below such that they already could not scatter the surface photons resonantly. 
Nonresonant scattering is strongly suppressed (Beloborodov & Thompson 2007) therefore 
these pairs freely flow along the magnetic field lines until they hit the surface of the neutron 
star where the energy is released. The observed hard X-ray persistent emission definitely 
comes from an optically thin, tenuous plasma; this suggests that the energy is released at a 
small depth, < 1 g/cm 2 . 

It was assumed by Thompson & Beloborodov (2005) and by Beloborodov & Thomp- 
son (2007) that the two-stream instability develops in a very narrow upper layer of the 
atmosphere so that collisionless processes are responsible for deceleration of the beam and 
plasma heating. The observed hard radiation was attributed to bremsstrahlung from this 
atmosphere. However, the two-stream instability is quenched when the beam is spread into 
a plateau-like distribution so that only one half of the initial beam energy is released in the 
optically thin domain. The beam is in any case stopped finally at a larger depth by collisions 
with the background particles. So only one half of the beam energy could be radiated in 
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the hard band. Of this, one half is directed towards the star where it is mostly absorbed 
and reradiated back in the soft band. One concludes that collisionless dissipation in the 
atmosphere converts only 1/4 of the total energy to the hard radiation and could not ac- 
count for the fact that most of the magnetars luminosity in the quiescent state is observed 
at e > 100 keV. Moreover, as demonstrated below, the inhomogeneity of the neutron star at- 
mosphere stabilizes the two-stream instability so that collisionless heating of the atmosphere 
is quenched. 

An alternative explanation is that the nonthermal tail is formed via resonant scattering 
of the thermal surface radiation at high altitudes, about 10 stellar radii, where the Landau 
energy is already nonrelativistic. This would mean that the persistent emission of the mag- 
netar is associated only with a small fraction of the magnetic field lines (those rising into 
the resonant region) whereas the currents flowing along the most of magnetic field lines and 
carrying most of the energy do not show up. Moreover, careful calculations show that the 
resonant scattering cannot reproduce the observed rising energy spectra of the persistent 
emission (Fernandes & Thompson 2007). 

In this paper, we reanalyze interaction of the fast particle beam with the surface lay- 
ers of the neutron star and show that, because of the extreme physics introduced by the 
ultrastrong magnetic field, a hot atmosphere could be formed via collisional processes and 
these difficulties avoided. It will be demonstrated below that even though these electrons 
appear at a relatively large depth, they could escape upwards forming a hot layer with the 
temperature T £_b/2 rs 1 2 MeV at the top of the cold atmosphere. The reason is 
that the Coulomb cross section sharply decreases when the electron energy becomes less 
than eb and the electron is restricted to move only along the magnetic field like a bead 
on a wire. These electrons could carry away a significant fraction of their initial energy if 
the atmosphere is fully ionized so that there are no ionization losses. We show that helium 
atmosphere satisfies these conditions. Such an atmosphere is expected in view of spallation 
of heavy nuclei bombarded by high-energy electrons and positrons. 

Bremsstrahlung radiation from this hot layer is hard enough to populate the magneto- 
sphere by electron-positron pairs. This pair corona is also heated via collisionless relaxation 
of the primary beam because the corona, unlike the atmosphere, is extended enough for the 
two-stream instability to develop. The observed hard X-rays are radiated via Comptoniza- 
tion in the corona and bremsstrahlung in the hot atmosphere. Our model predicts a hard 
radiation spectrum extending to the MeV band. In the band < 100 keV, Comptonization 
forms a power-law spectrum with the photon slope 1 < a < 2 as is observed. Because the 
corona is extended, most of radiation is not intercepted by the surface of the neutron star but 
rather escapes to the infinity. Therefore the observed persistent radiation from magnetars is 
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dominated by hard emission. 

The paper is organized as follows. In sect. 2, we analyze the two stream instability in the 
neutron star atmosphere and show that the instability is suppressed by the inhomogeneity 
of the atmosphere. In sect. 3, we consider deceleration of the fast particle beam via resonant 
Coulomb scattering and development of the electron-positron cascade in the atmosphere 
(Kotov & Kelner 1985, Beloborodov & Thompson 2007); we show that the initial energy 
of the beam is eventually converted into electrons with the energy something below eb 
at the depth of roughly 100 g/cm 2 . In sect. 4, we demonstrate that these electrons could 
escape upwards delivering most of the original energy into a hot layer at the top of the 
cold atmosphere. In sect. 5, we argue that hard radiation from this hot layer populates 
the magnetosphere with electron-positron pairs and this pair corona emits hard radiation 
extended into the MeV band. Conclusions are presented in Sect. 6. In Appendix A, we 
show that the energy the electron loses in the atmosphere for recoil of ions is small. In 
Appendix B, we demonstrate that helium is fully ionized in the atmospheres of magnetars. 
In Appendix C, we show that Comptonization of soft photons on mildly relativistic electrons 
in a super-strong magnetic field results in a power-law radiation spectrum with the photon 
index a > 1, like in the non-magnetized case. 



2. Inefficiency of the collisionless interaction of the beam with the atmosphere 

Thompson & Beloborodov (2005) and Beloborodov & Thompson (2007) suggested that 
the energy of the magnetospheric currents is dissipated when the downward electron-positron 
beam with the Lorentz factor is decelerated in a thin surface layer of the star where 
two-stream instability excites Langmuir turbulence thus providing the necessary relaxation 
mechanism. In this section, we check conditions for the development of the two-stream 
instability. One can conveniently express the plasma density in the atmosphere via the 
Thompson optical depth, r, as 

< 3 > 

where 

H= (l + Z)T = 26 (l + Z)T 7cm 
Am p g Ag u _ 5 

is the hydrostatic height of the atmosphere, g = 10 14 ' 5 gi4.5 cm-s -2 the surface gravity, A and 
Z the atomic weight and charge, correspondingly. The beam density may be expressed via 
the total energy release in the magnetosphere, L = 10 36 L36 erg/s, as 

n b = , J 3 =4xl0 15 ^ cm' 3 ; (5) 
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where 73 = 7/IO 3 . 

The basic regimes of the two-stream instability for a relativistic beam propagating in 
a non-relativistic plasma were studied by Fainberg, Shapiro & Shevchenko (1970). For a 
monochromatic beam moving with the Lorentz factor 7 along the strong magnetic field, the 
growth rate is written in the hydrodynamic regime, 

fn b \ l /zu p 

Khydr ~ W3 [-) "J (6) 

where 



«V = \— P) 

is the plasma frequency. This formula is valid provided all the particles in the beam are in 
resonance with the excited wave, 5v/c < k/uj p . For the beam with a large energy spread, 
5-y ~ 7, this condition is written as 

eyvi. (s) 

Making use of Eqs.([3]l6]), one can write this condition as 

7> 400 L AT ?}t\ V 72 - 0) 



1 + Z)T 7 L 



36/ 

In the opposed limit, only some fraction of the beam particles are in resonance with the 
excited wave (kinetic regime). Then the growth rate is 

Kkin ~ tt^ ~ u p . (10) 

n 7 d of 1 n 

Making use of these estimates, one easily finds that kH/c is very large in both regimes 
therefore at first glance, strong Langmuir turbulence should be excited. However, careful 
consideration shows that this is not the case because at k/uj p 1, the instability is easily 
stabilized by plasma inhomogeneity (e.g. Breizman & Ryutov 1974; Breizman 1990). Let 
us consider how inhomogeneity of the neutron star atmosphere affects development of the 
two- stream instability. 

The instability excites the Langmuir waves with the dispersion relation 

/ 3kT k 2 c 2 \ .... 
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where T is the plasma temperature, k the wave vector. Within the atmosphere with the 
characteristic height (j3J), the plasma frequency varies with the depth, z, as 

Su p _ 6z 

oo p - 2H- 1 ] 

The frequency of the propagating wave remains constant whereas the wave vector varies to 
satisfy the dispersion equation ffTTj) : 

IkT 

S U + - — 5k = 0. (13) 

m e c 

Here we take into account that for the resonance wave, uj/k « u p /k ~ c. The waves are 
amplified only if their phase velocity is close to the beam velocity, 

\oj - vk\ < k. (14) 

In the hydrodynamic regime, when the beam velocity spread is less than the width of the 
resonance, the condition (fl4"|) yields 5k < Khydr/c- Then one finds that the amplification 
stops after the wave propagates the distance 

5z = ^L^y± H . (15) 

m e c z Up 

A significant fraction of the beam energy could be dissipated provided 

KhydrOZ 



V 9 



> A; (16) 



where A = 10Ai is the logarithm of the ratio of the beam energy to the initial energy of 
the Langmuir oscillation (it is of the order of the Coulomb logarithm) and it is taken into 
account that the wave propagates with the group velocity 

du) 3kT .„ „. 

= Tk = ^ (17) 

Now one finds making use of Eqs.Q, (|T5l) and (|T7|) 

>^SZ = 13 x 1Q - 6 (l + ^7 /6 ^\ (lg) 

v e A r 1/6 Ai 
So the condition ffl6|) is violated by a large margin. 

In the kinetic regime, the condition (fill) reduces to 5k /k < 5v/c ~ (cy 2 )~ l , which yields 

5z = ~^H. (19) 



Then one finds 



v g A 



4.3 x 10 
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(1 + Z)T, 
Argils 



(20) 



Again, the condition f|T6l) is violated by a large margin. Therefore the two-stream instability 
does not develop in the atmosphere of the neutron star and the atmosphere could not be 
heated via collisionless dissipation. 

It will be shown below that a hot atmosphere could arise in the case under consideration 
but as a result of collisional interaction of the beam with the surface layers of the star. Then 
the characteristic scale of the density variation significantly increases and moreover, the hot 
atmosphere emits hard X-rays therefore the Lorentz factor of the beam decreases according 
to Eq.©. Substituting in Eq. ffTB"]) the characteristic temperature of the hot atmosphere, 
T = 10 10 K, and the Lorentz factor 7 = 100, one can see that the two-stream instability 
could develop in the hot atmosphere. Therefore eventually some fraction of the beam energy 
could be released via collisionless dissipation. However, at least one half of the energy is in 
any case dissipated and delivered into the hot atmosphere by collisional processes; it is these 
processes that determine the temperature of the hot atmosphere. 



Consider deceleration of a relativistic electron-positron beam within the surface layers of 
the magnetar. Let a relativistic electron move in plasma along the strong magnetic field. If 
the electron energy exceeds the Landau energy ([T]) , it efficiently scatters off ions into the first 
Landau level and then immediately falls back emitting a resonance photon; such a resonance 
Coulomb scattering is the main mechanism of deceleration of a relativistic electron-positron 
beam in the case under consideration (Kotov & Kelner 1985, Beloborodov & Thompson 
2007). The cross section for the transition of an electron with the Lorentz factor 7 and the 
momentum p from the ground to the j-th Landau level by scattering off an ion with the 
charge Ze is found by Bussard (1980) and Langer (1981): 



3. 



Collisional deceleration of the beam 
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Here Ot is the Thomson cross section, s' the final spin of electron and summation is over 
the final electron momentums p± = ±a/p 2 — 2m e c 2 jB/B q . Even high energy electrons jump 
predominantly on the first Landau level therefore one can neglect excitations of higher levels. 
The cross section for the transition — > 1 is plotted in Fig.(l). In the limit mc 2 7 ^> Eb, it 
is reduced to 

<no = ^Vr In (o.413^). (22) 

The electron with the Lorentz factor 7 > ^2B / B q jumps on the first Landau level 
having the Lorentz factor 71 = m e c 2r )/eB = 0.15 , ~fB l ^ 2 . The electron immediately falls 
back onto the background Landau level emitting a photon with the energy (0.5 -7- 1)eb 
as measured in the guiding center frame of the excited electron (frame moving with the 
Lorentz factor 7^; due to recoil, the energy of the emitted photon depends on the emission 
angle. After the deexcitation, the Lorentz factor of the electron remains on average the same 
because in the guiding center frame, photons are emitted forward and backward with equal 
probability. Thus the electron retains only a fraction £ = 71/7 ~ 0.15-B^ 1 / 2 of the total 
energy, most of the energy being taken away by a photon. 

The fate of the photon depends on its polarization. Two polarization modes could prop- 
agate in the magnetized vacuum; the so called ordinary mode is polarized in the plane set 
by the propagation direction and the background magnetic field whereas the extraordinary 
mode is polarized perpendicularly to this plane. If an ordinary photon is emitted, it is im- 
mediately converted into an electron-positron pair provided its energy exceeds the threshold, 
E > 2m e c 2 / sin 9. In magnetar's field, this condition is satisfied for most 9. In the frame of 
the scattered electron, the produced electron and positron have on average equal energies 
and move in opposite directions. In the laboratory frame, most of the energy is taken by 
the particle moving forward; only a fraction I/7 2 of the total energy is taken by the second 
particle. This means that the energy of electrons in the beam decreases only by a fraction £ 
in one scattering. If this were the only process of the beam-plasma interaction, the energy 
of the particles in the beam would decrease exp (r£ ) times after r scattering and the total 
number of scattering before the particle energy becomes less than the Landau energy would 
be large, r £ _1 ln^/eg) ~ 20 -f- 30. The full length of the avalanche, I, could be estimated 
by summing the free path lengths: 

l _ 1 AB r dn 

~ ~^ni ~ SBgZ^rrii J ln(0.4 7 2 " 2 £n 



23 h, 



3£B q Z 2 (jTn 
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Here rii is the number density of ions. For 7 ~ 100 -j- 1000, the corresponding Thomson 
depth is 



An extraordinary photon does not produce pairs directly; it first splits into ordinary 
photons and only then pairs may be produced. Then the avalanche proceeds further. As 
the energy per particle decreases in this case roughly twice in each step, the avalanche 
penetrates the depth significantly less than that of Eq. (1231) provided each emitted O-photon 
is converted into a pair. However, the energy of the photon, (0.5 1)£_b ~ 2-^3 MeV, is now 
distributed between two photons, both or one of them could fall below the threshold for the 
pair production. Then these photons either are converted into pairs in the Coulomb field of 
the nucleus or experience a Compton scattering off a background electron. 

The cross section for the pair production at a nucleus is 



where a is the fine structure constant. The Compton scattering of a photon moving along 
the superstrong magnetic field was studied by Gonthier et al. (2000). In this case, the 
resonance occurs only at the cyclotron fundamental e = (B/B q )m e c 2 = I2.B15 MeV. Above 
the resonance, the scattering cross section is close to the Klein-Nishina one 



The pair production at nuclei dominates at energies e > 250Z~ 1 m e c 2 . Note that this process 
comes into play only when the photon could not be converted directly into a pair. This occurs 
predominantly when E-mode resonant photons are emitted because these photons split into 
O-photons with lesser energies so that these O-photons could fall below the threshold for 
the direct pair production in the magnetic field. In this case, two pairs are produced in 
each step therefore the energy per particle decreases roughly four times in each step. This 
means that pairs are produced at nuclei only in the first few generations; after this, Compton 
scattering dominates. The corresponding depth varies from rz ~ 200/Z at small Z, when 
the transition energy is high and the Compton scattering comes into play already after two 
generations, to Tz ~ 800/Z, when Z is large and three or four generations are necessary. 

At the Compton stage, the avalanche proceeds further. The distribution of the scattered 
photons in their energies, e', is 




(23) 




(24) 




(25) 



3 m e c 2 de' e e' 




) 



2 




dac = 7:<?t 



- 2 



(26) 
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where 

e , 
1 + 2e/m e c 2 ~ £ ~ £ - 

The scattered photon takes on average a fraction 4/[3 ln(2e/m e c 2 )] ~ 0.2 of the initial en- 
ergy. The photon is directed at the angle ~ ^/4m e c 2 /e to the magnetic field therefore it 
is immediately converted into a pair. The rest 0.8 fraction of the energy is taken by the 
recoil electron, which emits a resonant photon via Coulomb scattering off an ion. If this 
photon is in O-mode, it is immediately converted into a pair; if it is in E-mode, it splits 
producing two photons, the energy of each of them being roughly (1/2) • 0.8(1 — £) ~ 1/3 of 
the initial energy. These two photons either produce pairs directly or, if they are below the 
threshold for the direct pair production, experience Compton scatterings and the process 
repeats again. The Compton scattering comes into play predominantly if an E-mode photon 
was emitted because a resonant O-photon is converted into a pair immediately whereas the 
E-photon splits into two O-photons, which could fall below the threshold of the direct pair 
production. Therefore in the Compton channel, the energy per particle decreases by the fac- 
tor of 3 in each generation. Then the energy per particle becomes less than Eb in a few steps 
and the avalanche stops. As the Compton cross section grows rapidly with decreasing of the 
photon energy and cyclotron resonances only increase the cross section, the full length of the 
avalanche is determined by the free path of a photon with the initial energy e = 250Z~ 1 m e c 2 . 
The corresponding depth varies from tq ~ 100 at Z = 1 to t c ~ 10 at Z = 26. We will argue 
in the next section that for our model to be self-consistent, the avalanche should develop in 
the medium composed from light elements. Spallation of nuclei by relativistic electrons and 
positrons from the avalanche seems to result in formation of a helium atmosphere; below we 
will consider only this case. Then the total depth the avalanche penetrates is estimated as 
T = Tz + t c ~ 300/Z. 

In real avalanches, ordinary and extraordinary photons are emitted alternately with 
approximately the same probability. As the Coulomb cross section ( 12TI) is significantly larger 
than (1231) and (124")) . the longest parts of the chain is associated with emission of extraordinary 
photon with subsequent photon splitting. On the other hand, the total number of generations 
is also determined by emission of extraordinary photons because the energy per electron 
decreases significantly when the photon splits. Therefore the overall length of the avalanche 
is roughly the same as if only extraordinary photons were emitted. So finally the energy of 
the primary beam is converted into electron-positron pairs with the energy ~£g/2~l-h2 
MeV at the depth of 

r = arZriil ~ 300/Z. (27) 
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Fig. 1. — Cross-section of the Coulomb excitation of the first Landau level; asymptotics ( 1221) 
is shown by dashed lines. 
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Fig. 2. — Coulomb cross-section for the electron on the ground Landau level, B = 10 15 G. 
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4. Formation of a hot atmosphere 



Now let us consider the fate of the newly formed pairs with the energy less than eb- Their 
motion is purely one-dimensional and conservation of the energy and momentum implies 
that two colliding particles of equal mass may only exchange their energy and momenta. 
Therefore electron-electron collisions do not change the state of the system and may be 
neglected. Collision of a "hot" positron with a "cold" background electron results in a 
"cold" positron and "hot" electron. Collision frequency of cold particles is high therefore 
the cold positron does not escape but rather eventually annihilates with some background 
electron. However, the hot electron takes most of the energy in this case therefore most 
of the beam energy is eventually stored in electrons with the characteristic energy ~ 1 4- 2 
MeV. They diffuse within the medium colliding with the background ions. 

Collision of a one-dimensional electron with the ion may result either to the forward 
scattering, after which the electron energy remains the same to within a small recoil, or to 
reflection from the ion. The cross section for the Coulomb reflection is found from Eq. (|21|) 

as 



It was shown in the previous section that the avalanche penetrates the Thomson depth fl27j) 
producing electrons with the energy ~ £b/2. One sees that for these electrons, a^riil = 
0.3B^ 2 < 1 so they could easily escape upwards taking away most of the energy of the 
avalanche. The cross- section of the forward scattering is larger than f[2"9~l) but because recoil 
is small, one can neglect this process (see Appendix A). 

The above consideration assumes that the escaping electron does not lose energy on 
ionization so that the plasma in the atmosphere is fully ionized. It is shown in the Appendix 
B that under the condition of interest, helium is fully ionized. One can assume that the 
layer the avalanche penetrates composed mostly of helium because the avalanche electrons 
and gamma-quanta destroy heavy nuclei (Beloborodov & Thompson 2007). The helium is 
presumably preferable to hydrogen because does not require transformation of neutrons into 
protons. 





this cross section is plotted in Fig. 2. For p = 0.5eb/c = a/ B/2B q 



m e c one gets 



(T 00 = 0.6Z 2 a T (B q /B) 2 = 0.001Z 2 a T B^ 2 . 



(29) 



One should check how much energy the electrons lose on bremsstrahlung before they 
escape. In the non-magnetized medium, a relativistic electron loses energy on bremsstrahlung 
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after passing the depth, in Thomson units, 

T ° br = 3aZ(ln2 7 -l/3)' (30) 

Unfortunately bremsstrahlung in the superstrong magnetic field has not been calculated yet 
however simple quasiclassical estimate shows that the rate of the energy losses decreases in 
this case. 

In classical electrodynamics, the energy radiated by a relativistic electron may be written 
as (e.g. Landau & Lifshitz 1995) 

58 = % f 76[(v ' w)2 + cW /7 2 ]^; ( 31 ) 

where w is the electron acceleration. From the equation of motion 

m e — : — = eE 
dt 

one finds Wn = eE\\/m e c~f 3 and w± = eE±/m e cj where the subscripts || and J_ refer to 
the projections onto the direction of velocity and onto the perpendicular direction. In the 
non-magnetized case, both components of the acceleration are presented however the first 
term in the square brackets in Eq. (l3T|) is 7 2 less than the second one and may be neglected. 
For the electron passing the nucleus at the distance r one gets assuming that the motion is 
straightforward 

vrZ 2 e 4 7 2 

68 = Ylr*mW (32) 



Integrating the obtained relation over r from r m j„ = firy /m e c defined from the condition that 
the energy of the emitted quanta becomes equal to the electron energy, one gets Eq. (l30l) to 
within a numerical factor and the logarithmic term. 

In the superstrong magnetic field, the electrons move only along the field therefore only 
the first term in the square brackets in Eq. fl3T|) remains. Simple calculation shows that in this 
case 68 is 37 2 less than that of Eq. fl32|) . One can expect that the full bremsstrahlung cross 
section in the superstrong magnetic field decreases by the same factor. Of course quantum 
treatment of the process is necessary to rigorously justify this conclusion, but assuming that 
the above consideration is basically correct, the electron in the superstrong magnetic field 
loses the energy on bremsstrahlung after passing the depth r br = 37 2 TQ r . For electrons with 
the energy ~ £b/2, one gets r br = 6000Z _1 5i5, which is significantly larger than the depth 
( 1271) the avalanche penetrates. This means that one can neglect bremsstrahlung energy 
losses. 
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So hot electrons freely escape upwards taking away most of the beam energy. Charge 
neutrality implies that the necessary amount of background ions rises together with the 
electrons so that a hot atmosphere with the temperature T^£#/2~l-h2 MeV and the 
characteristic height H = 26(1 + Z)T W /Agi^ m is formed. 

5. Pair corona 

The hot atmosphere is cooled via bremsstrahlung; radiation is dominated by photons 
with the energy ~ T ~ 1 4- 2 MeV. One half of this radiation is directed towards the star; it 
is absorbed and eventually is reradiated as thermal emission. The other half of radiation is 
radiated away. In the magnetosphere, hard photons could be converted into electron-positron 
pairs. The pair production rate may be written as N = (L/m e c 2 where the numerical factor 
£ < 1/2 takes into account uncertainties in the radiation spectrum and the field geometry. 
Within the light travel time, R*/c, the magnetosphere will be filled by pairs with the density 
n = 8 x 10 18 C-^36 cm -3 . The characteristic inhomogeneity scale of the corona is of the 
order of the star radius, which is large enough for the two-stream instability to develop in 
the corona. The condition for the instability f|T6|) is satisfied provided the left-hand side of 
Eq.ffTBl exceeds unity. Substituting in this relation H by and making use of Eqs.flH]), ((7j) 
and the above estimate for the density of the corona, one gets 

A T 1 / 2 
^kinOZ _ ^q3 -^36 

v g A ~ C 1/6 7 3 8/3 Ai' 

Note that the condition (jSJ) is satisfied in the corona at any reasonable parameters therefore 
the instability is hydrodynamic. So after the corona is formed, the primary beam would 
experience collisionless relaxation. Then about one half of the beam energy is spent on 
heating of the corona whereas the other half will be delivered to the surface of the star where 
the the electron-positron avalanche is developed and the energy of the beam is eventually 
transferred to electrons with the energy ~ £b/2 ~ 14-2 MeV each. It was argued in the 
preovious sections, that these electrons escape upwards delivering most of their energy into 
the tenious, hot atmosphere. 

The corona is efficiently cooled by Comptonization of the thermal emission from the 
surface. As the energy of the photons is well below the Landau energy, only O-mode radiation 
is scattered (polarized in the plane set by the direction of propagation and the magnetic field). 
Thermal energy stored in the star interior is transferred to the surface and radiated away 
by E-mode photons because their opacities are (Bmc 2 / B q e) 2 » 1 times less than those for 
O-mode photons. Soft O-mode photons could be emitted only if surface layers of the star 
are heated. As hard radiation from the hot atmosphere illuminates the underlaying cold 
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atmosphere, some fraction of this radiation will be absorbed and reradiated in the soft band. 
It is this thermal O-mode radiation that could be a soft photon source for Comptonization. 

The scattering cross section for the O-mode photons is a = a? sin 2 9' where 9' is the 
angle between the propagation direction and the magnetic field in the proper electron frame. 
The relativistic electron "sees" radiation at the angle 9' ~ I/7 therefore the cooling rate 
decreases 7 2 times as compared with the non-magnetized case and may be estimated as 
(Beloborodov & Thompson 2007) 

(f ) c = < 33 > 

where £ is the energy of the electron (assumed to be larger than m e c 2 ), U the radiation 
energy density. Writing the radiation energy density as U — L/AnRlc, one finds that the 
electron cooling time 

t = Ax lO-^McV^ae 1 (34) 
is comparable with the light travel time. 

The observed radiation is a superposition of the bremsstrahlung radiation from the 
hot atmosphere and the Comptonization radiation from the corona; the radiation spectrum 
extends to the characteristic particle energy of l-i-2 MeV. Bremsstrahlung has a flat intensity 
spectrum (the photon index 1) at e <C T whereas unsaturated Comptonization produces a 
power-law spectrum with the photon spectral index a > 1 depending on the parameters of 
the system (see Appendix C). As the luminosity of the corona is larger than the thermal 
luminosity of the surface, which provides soft photons for Comptonization, the slope should 
be hard enough a < 2. Therefore the low- frequency part of the spectrum is dominated by 
Comptonization and could exhibit a variety of spectral indices in the range 1 < a < 2 as 
is observed. The observed high pulsed fraction is naturally explained by the fact that the 
scattering cross section for the ordinary mode photons is highly anisotropic. 

When the pairs fill the magnetosphere, the energy release stops because the pair plasma 
shorts out the induction electric field. Energetic primary particles, which have already filled 
the magnetosphere, disappear at the star's surface for about the light travel time. At the 
next stage, the hot atmosphere and the corona are cooled. The cooling time is of the order of 
the light travel time. An important point that positrons could not survive for the larger time 
because they are annihilated when hitting the surface of the star. When a positron falls from 
the corona onto the surface, it exchanges energy with an electron in the cold atmosphere 
and annihilate^]. 



Note that even though the annihilating particles are cold, the annihilation line is not formed. In the 
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The corona is expected to be extremely unsteady. Even though most of the energy 
is taken away by the electron, which reflects from an ion and goes upwards, this electron 
cannot rise directly into the corona because of the charge neutrality and would remain in 
the hot atmosphere. The corona is replenished by the hydrodynamic expansion of this 
hot atmosphere, and cooled by Compton losses. The Compton losses may or may not be 
compensated by collisionless interaction with the primary beam. Therefore the corona may 
be depleted by Compton cooling; in any case, we expect that half the primary beam energy 
makes it down to the hot atmosphere, so that new matter is constantly being added to the 
corona, and this very likely leads to a non-steady situation. Because the coronal matter can 
propagate the currents, the displacement currents and attendant primary beam are switched 
off, while Compton cooling continues. When the temperature of the hot atmosphere falls 
below m e c 2 , the pair production stops and the plasma density in the corona falls below the 
critical value necessary to maintain the magnetospheric currents. Then the displacement 
current arises again and the next cycle of energy release starts. So one can expect strong 
fluctuations of the radiation at the timescale of At > 10~ 4 s. The possibility exists of 
revealing them by analysis of the photon statistics. If a source emits radiation in separate 
bursts of the characteristic duration At = 10~ 3 At_ 3 , the probability distribution for a pair 
of photons to be detected within the time interval t should differ significantly from a Poisson 
distribution at t ~ At. For a collecting area of 10 3 A 3 cm 2 , and a photon flux of 10 _3 .F_ 3 cm~ 2 , 
an observation interval of 10 5 tsS should contain 10 5 A^F^t^ photons and, for Poisson arrival 
statistics, 10 2 A 3 F 2 3 t 5 At _ 3 pairs of photons arriving within At of each other. The non-steady 
nature of the hard coronal emission (this emission in fact dominates the thermal surface 
emission already at a few keV), which produces deviations from Poisson arrival statistics, is 
therefore detectable with a sufficiently powerful detector and large exposure times. 

6. Conclusions 

We considered dissipation of the energy released in the magnetosphere of the magnetar 
in the course of slow relaxation of non-potential magnetic fields. A basic physical picture was 
proposed by Thompson et al. (2002) and recently elaborated by Thompson & Beloborodov 
(2005) and Beloborodov & Thompson (2007). When the plasma density in the magneto- 
sphere falls below a critical value necessary to maintain the magnetospheric currents, an 
induction electric field arises and initiates an electron-positron avalanche resembling that 



superstrong magnetic field, only longitudinal (along the field) component of the momentum is conserved 
therefore the two-photon annihilation results in photons with generally different energies, the annihilation 
spectrum being extended from to 1 MeV (Kaminker, Pavlov & Mamradze 1987). 
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in pulsars. The fast electron-positron flow hits the surface of the star where the released 
energy is dissipated. The observed very hard spectra of the persistent X-ray emission from 
SGRs and AXPs (Kuiper et al. 2004, 2006; Mereghetti et al. 2005, Molkov et al. 2005) 
imply that most of the energy is released in a very hot and tenuous plasma. It was assumed 
by Thompson & Beloborodov (2005) and Beloborodov & Thompson (2007) that the flow 
loses a significant fraction of its energy in a thin surface layer where the two-stream insta- 
bility develops so that the plasma is heated by collisionless processes. Collisionless heating 
is balanced by bremsstrahlung radiation and the equilibrium temperature about 100 keV is 
achieved. 

We reanalyse interaction of the fast plasma flow with the surface of the magnetar and 
conclude that this mechanism is incapable of heating the atmosphere because strong density 
gradient in the atmosphere of the neutron star suppresses the two-stream instability. We 
propose, rather, that a hot, tenuous atmosphere/corona could arise due to specific properties 
of Coulomb scattering in the superstrong magnetic field, mainly due to one-dimensional 
character of the electron motion. Within the upper layers of the neutron star, the flow energy 
is transferred, via an electron-positron avalanche, to electrons with the energy something less 
than the Landau energy (~ 1 -f- 2 MeV in the magnetar's magnetic field). This occurs at 
a significant depth, ~ 100 g/cm 2 , however these electrons are not thermalized but rather 
escape upwards taking away most of the initial flow energy. The reason is that collisions 
between electrons in one-dimension do not result in relaxation; after the collision, the two 
energies of the two electrons are the same as the initial energies. On the other hand, collisions 
with ions result only in nearly elastic reflection. These electrons form a hot atmosphere (the 
necessary amount of ions accompany the electrons in order to maintain charge neutrality) 
with the temperature ~ l-=-2 MeV, which is an order of magnitude larger than in the model by 
Thomson & Beloborodov (2005) and Beloborodov & Thompson (2007). Hard radiation from 
this atmosphere is generated via bremsstrahlung. Pairs are easily produced in this radiation 
field; they fill the whole magnetosphere forming a hot corona. Collisionless interaction of 
the primary beam with the pair plasma in the corona heats the pairs even more; they are 
cooled by Comptonization so that the overall spectrum of the source is a superposition of 
the bremsstrahlung radiation from the hot atmosphere and a Comptonization radiation from 
the corona. 

The extended corona radiates in all directions so that only a small, < 1/2, fraction of the 
radiated energy is intercepted by the surface of the star. Therefore the observed luminosity 
is dominated by the hard radiation. The spectrum is extended to MeV band. Unsaturated 
Comptonization generates a power-law spectrum, which is generally steeper than the flat 
bremsstrahlung spectrum therefore radiation from the corona dominates in the range < 100 
keV. The photon spectral slope is 1 < a < 2, as is observed. The energy release in the 
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magnetosphere occurs spasmodically at the time scale of at least the light travel time: the 
pairs short out the induction electric field in the corona and the energy release stops until 
the corona cools down, then the charge starvation necessitates the displacement current and 
the process starts again. Therefore one can expect strong fluctuations of the radiation at 
the time scale of > 10~ 4 s. 

We are grateful to Rashid Shaisultanov for help. Y.L acknowledges support from the 
German-Israeli Foundation for Scientific Research and Development. D.E. acknowledges 
support from the United States - Israel Binational Science Foundation and from an Israel 
Science Foundation Center of Excellence Grant. 



Appendix A. Electron-ion collisions: recoil effect 



When the electron with the energy less than e b passes an ion, the energy is transferred 
only due to recoil effect. In magnetar's field, most of the ions are in the ground Landau state 
and a scattering occurs if the ion makes a transition to the first level with the energy 

ZeB Z 

e B i = = 5— .Bis keV; 

rriiC A 

where rrii = Am p is the ion mass. The cross section for collisional transitions between the 
ion Landau levels was found by Langer (1981); for the transition between the ground and 
the first levels it is written as 



a^Za T mlc y {\nK- 0.577); 
pp' 



8B 



where 



A - 1 = | 



^Y-(7-Y) 2 

m e c ) 



The conservation laws imply 

m e c7 + rriiC = m e cy' + mjc 2 + (p - p') 2 + 2e B imi\ 

which yields 



and 



7 - 7 = ^Bi\ 

£Bi ( ra e C7 _ I 
rriiC 2 VP 4 
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As the electron looses only a small fraction of the energy in a scattering, 7 - 7' <C 7, one 
can conveniently define the effective cross-section for the energy loss 

Z 2 

a = [(7 -i)h\a = 2 x 1(T 3 — [1 + 0.08 ]n(A 2 /ZB 1B )]a T . 

It was shown in Sect. 3 that the electron-positron avalanche stops at the depth ( 1271) forming 
eventually a cloud of electrons with the energy ~ £b/2. The electron with this energy loses 
the fraction of its energy 

A T? 1 

= an t l = 0.1 m [\ + 0.081n(A 2 /^ 15 )] 

when rising from the depth ( 1271) . This estimate assumes that the electron moves straight- 
forwardly but not diffuses upwards. This is justified because the layer of the depth f[2"Tj) is 
transparent for the Coulomb reflection (see the estimate (129)) ). So the fraction of the energy 
the electron loses for recoil process is small. 



Appendix B. Ionization equilibrium. 

Ionization equilibrium in the super-strong magnetic field is still a subject of intense 
research, see recent reviews by Lai (2001) and Harding & Lai (2006) and references therein. 
Here we present rough estimates for the helium plasma. The ionization energy of the 
hydrogen- like ion is 



Q = 0.16Z 2 



In. 



h 3 B ^ 



a.u. 



m 2 e 3 cZ 2 

For the helium ion, one gets Q = 2.39(1 + 0.086 In B^) keV. When the atom moves, the 
ionization energy decreases so that the above value gives the estimate of the ionization 
temperature from above. 

The ionization equilibrium He ++ <=^He + + e is written as 

n e n ++ _ Z e Z ++ ^ 

n + Z + 

where n e , n ++ and n + are the number densities of free electrons, He ++ and He + , correspond- 
ingly, Z e) Z ++ , and Z + their partition functions. The partition function of the strongly 
magnetized electrons is 
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The ratio of partition functions of He ++ and He + is dominated by the ionization energy and 
may be presented as 



exp - 



Q 
T 



Now the temperature of ionization (when n ++ = n+) is found as 

3<=R f m..T\ 1/z 1 

T in „ = Q < In 



2imhc 



( m e T \ 
\2irh 2 J 



1 7 x 1Q6 (l + 0.091n^ 15 ) 2 

1 + 0.06 hi (£15/^14.5) 



This means that helium is fully ionized in magnetar's atmosphere with the temperature 
T = 0.5 4- 1 keV. 



Appendix C. Comptonization in a superstrong magnetic field 

Here we demonstrate that Comptonization of soft photons in the superstrong magnetic 
field results in a power-law spectrum with the photon index a > 1, like in the non- magnetized 
case. Compton scattering of soft photons on hot electrons results in a photon flux over the 
spectrum from the initial energy, e , to the spectral region e ~ T. If the optical depth 
of the source is large enough, the photons are accumulated at e ~ T and the equilibrium 
Bose-Einstein spectrum, ./Vbe = {exp[(r^ + e)/T] — is formed. This regime is called 
saturated Comptonization. In the medium of the moderate optical depth, photons are not 
accumulated but rather escape and therefore a power law spectrum could be formed in the 
range Eq <^ e <^ T. Here we show that the same occurs also in the magnetic field so strong 
that the electrons populate only the ground Landau level. For nonrelativistic temperatures, 
Comptonization in the superstrong magnetic field was studied by Lyubarskii (1987a,b) in the 
Focker-Plank approximation. Here we allow relativistic temperatures but restrict ourselves 
only to low photon energies when one can neglect recoil. 

Let us first consider scattering on electrons moving along the magnetic field with some 
momentum p; the number density of this electrons is f(p)dp; where f{p) is the electron 
distribution function. The kinetic equation for photons is easily written in the proper frame 
of these electrons as 

( J7 + c1 '^) e', V) = cf'dp* J deWie', - e')^[n(r', s\, l'J - n(r', e', 1% 

where 1 is the direction of propagation of photons, n(r, e, 1) the phase density of photons and 
prime marks quantities measured in the proper electron frame. The scattering cross-section 

is ^3 

- = — a T sin 2 9' sin 2 9\ ; 
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where 9 and 9± is the angles between the photon direction and the magnetic field before 
and after the scattering, correspondingly. One can transform this equation to the laboratory 
frame taking into account that the distribution functions are relativistic invariants, f'(jJ) = 
f(p), n'(e', V) = n(e, 1), as well as the expressions ededVt and + lj^). Summation over 
all electrons yields 

din din 3 f 

sin 2 6* sin 2 9 1 f(p)dpde 1 dQ 1 



x5[e(c — vcos9) — e\(c — vcos9\)\- 



7 6 (c — v cos9)(c — v cos ^i) 3 ' 
This equation describes Comptonization of photons with the energy larger than the energy 
of the seed photons, e , but small enough for the recoil effect to be neglected, e <^ e <ti 
min(T, m e c 2 /T). 

In the steady state case, d/dt = 0, solution to this equation has a power law form 

n(r,e,l) = J(r, 1)^( 2+Q ); 
where the spatial function J satisfies the equation 

J(r,h)- J(r,l) 



/JM) 3 5 



c — V COS t 



f cos^ 



sin 2 9 sin 2 9\f {jp)dpdVL\ 



7 6 (c — f cos 6*)(c — f cos ^i) 4 

(Bl) 

The photon power index, a, could be found as an eigenvalue of this equation. It is determined 
by the electron distribution function and by the geometry and the optical depth of the 
source. It is beyond the scope of the present paper to solve this equation (solution of a 
similar problem for a nonmagnetized plasma is given by Titarchuk & Lyubarskij (1995)). 
Let us only demonstrate that one can expect a > 1. 

In the infinite homogeneous medium, the left-hand side of Eq.(Bl) is zero therefore one 
gets the equation 

/ - - <« - y:'^^ - * (-) 

where \x = cos9. This equation has an evident solution a = —2, J(fi) = const, which is 
nothing more than the low frequency part of the equilibrium Bose-Einstein spectrum. We 
are interested in a solution with a non-zero photon flux over the spectrum; in the infinite 
medium, such a solution implies that except of the soft photon source at e = Eq, there is 
a sink at some large enough energy e sink ; then Eq.(B2) describes the region e < e < £ sink . 
In a non- magnetized plasma, the solution with a non-zero photon flux over the spectrum is 
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n oc e -3 (Kats, Kontorovich & Kochanov 1978) so that the intensity spectrum is flat. In 
order to see that the same spectrum (a — 1 in our notations) satisfies also Eq.(B2) note that 
at a — 1, the integrand in the left-hand side of this equation is antisymmetric with respect to 
exchange /j «-> fa. Therefore the integral from the left-hand side of Eq.(B2) over /x vanishes 
identically. This means that a finite linear homogeneous set of equations, which could be 
obtained from Eq.(B2) by discrete approximation of the integral, is linearly dependent and 
therefore it has a nontrivial solution. ( The formal proof. Denote the linear operator in the 
left-hand side of Eq. (B2) at a = 1 as C and introduce the standard notation for the scalar 
product of functions (ip, <p) = j\ tpc/xifi. Now one can write that (e, CJ) = for e = const 
and an arbitrary J. Then (C*e, J) = so that e = const is a nontrivial solution to the 
conjugate equation C*e = 0. In this case, the equation CJ = also has by the Fredholm 
alternative a nontrivial solution.) 

Thus the spectrum with the slope unity is formed in the infinite medium; then all 
photons produced at e reach £ s i n k- In the case of a finite optical depth, the photons escape 
so that the spectral photon density should decrease with the frequency faster than in the 
infinite medium. Therefore unsaturated Comptonization in the super-strong magnetic field 
generates power law spectra with the slope a > 1 like in the non-magnetized plasma. 
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